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AiirHACT  m  «!##  ir  mw  ^ 

There  is  an  ongoing  interest  in  simulating  nonideal  blast  environ¬ 
ments  for  nuclear  effects  research.  In  particular,  one  would  like 
to  be  able  to  impose  peaked  blast  waves  on  real  ground  surfaces 
and  experimentally  measure  the  ensuing  dusty  airblast  environment. 
Proposed  here  is  a  gaseous  detonation  tube  blast  simulator.  A  dis¬ 
posable  (or  reusable)  shock  tube  would  be  constructed  on  an 
in-situ  ground  surface  of  interest.  The  tube  would  be  sealed  and 
filled  with  a  detonatable  gas  mixture.  When  a  planar  detonation 
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wave  is  initiated  at  one  end  of  the  tube,  it  induces  a  peaked 
blast  wave  which  expands  self-similarly  with  time — the  longer  the 
detonation  run  distance,  the  longer  the  blast  wave  duration. 
Similarity  analysis  of  such  a  wave  (which  consists  of  a  constant- 
velocity  Chapman-Jouguet  detonation  followed  by  an  adiabatic 
rarefaction  wave  expressed  in  terms  of  a  Riemann  characteristic) 
results  in  a  closed-form  analytic  solution  for  the  flow  field  time 
history.  It  is  shown  that  the  static  and  dynamic  pressure  wave^ 
forms  associated  with  this  detonation  give  a  high  fidelity  simula¬ 
tion  of  a  nuclear  surface  burst.JS^e  can  achieve  peak  overpres- 
>  sures  from  14  to  55  atm  and  peedc  dynamic  pressures  from  6  to  35 
atm,  depending  on  the  detona table  gas  selected.  Detonation  run 
lengths  from  25  to  100  m/KTsB^/^  dre  required  for  kiloton-level 
simulations.  The  principal  simulation  deficiency  involves  the 
detonation  temperature  (''^3000K)  which  is  considerably  larger  than 
its  airblast  counterpart  ('n<1000K)  .  Small  dust  particles  have  the 
potential  for  melting  and/or  vaporizing,  due  to  convective  heat 
transfer  from  the  detonation  products,  while  in  the  nuclear  case 
such  effects  are  caused  by  radiative  heating  from  the  fireball. 
These  effects  can  be  ameliorated  somewhat  by  choosing  air  as  the 
oxidizer,  and  by  employing  large  diameter  dust  particles. 
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1.  INTRODUCTION 


There  is  an  ongoing  interest  in  simulating  nonideal 
nuclear  blast  environments  for  weapons  effects  research. 
Gilette  et  al.  (1982)  have  used  a  portable  wind  tunnel  to 
measure  the  steady  flow  erosion  characteristics  of  desert 
soils.  We  would  like  to  do  similar  tests  with  peaked  blast 
waves  imposed  on  real  ground  surfaces  to  diagnose  experi* 
mentally  the  dusty  airblast  environment.  Dynamic  similitude 
arguments  applied  to  such  dusty  flows  lead  to  the  require¬ 
ment  for  large  explosive  yields. 

Solid  explosive  charges  have  been  used  to  create 
surface  burst  (and  more  recently,  height-of-burst)  blast 
environments  for  yields  as  large  as  a  kiloton  nuclear  equiva¬ 
lent.  Such  tests  are  expensive,  infrequent,  and  at  present, 
they  give  a  high'fidelity  simulation  for  peak  pressures  only 
below  about  80  psi.  Large  shock-tube  facilities  (such  as 
Sandia's  thunderpipe  and  the  DASACON  conical  tube)  are  usually 
limited  to  peak  pressures  of  tens  of  psi,  and  must  necessar¬ 
ily  use  disturbed  soil  seuoples.  The  dynamic  airblast  simula¬ 
tor  (DABS)  concept  could  be  used  for  creati.ng  kiloton-level 
blast  environments  (Renick,  1979)  with  peak  pressures  of  as 
high  as  600  psi,  but  test  times  are  short  (due  to  the  solid 
explosives  combustion  products)  and  such  tests  have  been  too 
expensive  for  parametric  studies.  Hence,  it  appears  that  the 
existing  blast  simulators  all  have  certain  limitations  for 
this  application.  There  is  a  need  for  a  simulator  which  is 
capable  of  producing  high-fidelity  nuclear  blast  environments 
with  peak  pressures  of  hundreds  of  psi  and  kiloton-level 
yields  on  in-situ  ground  surfaces.  Tests  run  with  such  a 
simulator  should  be  inexpensive  enough  that  parametric  studies 
can  be  performed. 
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A  new  concept  which  may  fulfill  such  requirements  is  the 
gaseous  detonation  tube  blast  simulator  shown  in  Figure  1. 

A  reusable  or  disposable  shock  tube  would  be  constructed 
over  an  in-si tu  ground  surface  of  interest.  The  tube  could, 
for  example,  be  constructed  of  concrete  sewer  pipe  of  a  semi¬ 
circular  cross  section,  with  a  soil  overburden  to  contain 
the  blast.  The  tube  Would  be  sealed  and  filled  with  a 
detona table  gas  mixture.  The  firing  of  the  detonator  on  the 
end  wall  would  initiate  a  detonation  wave  which  would  propa¬ 
gate  down  the  tube  at  a  constant  velocity,  known  as  the 
Chapman- Jouguet  (CJ)  velocity. 

The  detonation  wave  induces  a  peaked  blast  wave  flow 
field.  Figure  2  presents  the  flow  field  space  distribution 
of  a  planar,  constant  velocity  CJ  wave.  Pressures  (p) , 
densities  (p),  internal  energies  (e)  and  velocities  (u)  have 
been  nondimensionalized  by  their  values  just  behind  the 
detonation  wave  (denoted  by  subscript  n),  while  distances 
are  nondimensionalized  by  the  detonation  front  location. 

Such  a  wave  expands  self-similar ly;  that  is,  along  lines  of 
r/t  -  constant  the  flow  field  properties  remain  constant. 

Also  shown  in  Figure  2  is  the  flow  field  corresponding  to  a 
point  explosion  (PE) .  The  pressure  and  velocity  fields  for 
the  CJ  and  PE  solutions  are  qualitatively  similar,  while  the 
density  and  internal  energy  distributions  are  somewhat  dif¬ 
ferent  due  to  the  high- temperature  fireball  in  the  PE  case* 

The  remainder  of  the  paper  investigates  the  feasibility 
and  characteristics  of  such  a  simulator.  Section  describes 
the  similarity  analysis  for  one-dimensional  CJ  detonation 
waves,  and  gives  a  closed-form  analytic  expression  for  the 
planar  case.  Time  histories  for  the  planar  CJ  solution  are 
then  compared  and  matched  in  Section  XZX  to  a  1-KT  nuclear 
surface  burst  over  an  ideal  surface.  Conclusions  and  recom¬ 
mendations  are  offered  in  Section  XV. 
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x^ure  1.  Oetonafcion  tube  airblast  simulator 
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II.  SIMILARITY  ANALYSIS 


1.  GOVERNING  EQUATIONS 

Detonation  waves  are  nonsteady  gas  dynamic  flow  fields 
bounded  by  shock  waves  that  release  chemical  energy.  Among 
the  first  to  investigate  the  properties  of  such  waves  were 
Chapman  (1899)  and  Jouguet  (1905-6)  who  found  experimentally 
that  detonation  waves  propagate  at  a  constant  speed.  The 
fundamental  theory  of  detonation  waves  may  be  found  in  the 
books  by  Courant  and  Friedrichs  (1948)/  Stanyukovich  (1959), 
Zeldovich  and  Kompaneets  (1960) ,  emd  in  the  treatise  by  Taylor 
and  Tankin  (1958).  Analysis  of  TNT  explosions  is  contained 
in  the  work  of  Taylor  (1950).  Kuhl  and  Seizew  (1978)  com¬ 
pared  the  flow  field  of  self-similar,  planar,  cylindrical 
and  spherical  detonations  in  solid  and  gaseous  explosives  for 
the  constant  velocity  case,  while  Oppenheim  et  al.  (1972) 
parametrically  studied  all  possible  self-similar  blast  waves 
bounded  by  a  strong  CJ  detonation  front. 

Here  we  analyse  flow  fields  associated  with  planar 
detonations  under  the  following  idealized  conditions; 

•  The  ambient  medium  ahead  of  the  wave  is 
quiescent,  uniform,  and  unbounded. 

e  The  detonation  propagates  at  constant  velocity. 

e  The  flow  field  behind  the  front  is  adiabatic 
(no  afterburning)  with  a  constant  isentropic 
sound  speed  modulus. 

e  The  reaction  zone  is  negligibly  thin. 

•  The  detonate  vh  initiation  distance  is 
negligibly  small  compared  to  the  front 
radius. 

•  The  flow  field  is  plane  symmetrical  and 
unsupported  (e.g.,  by  a  piston). 


As  established  both  theoretically  and  experimentally, 
such  detonation  waves  propagate  at  a  constant  velocity  (the 
CJ  velocity) ,  and  satisfy  the  CJ  hypothesis  (the  wave  velocity 
equals  the  gas  velocity  plus  the  sound  speed) .  The  front 
trajectory  is  given  by 


r 


n 


w  t  . 
n  n 


(1) 


The  flow  field  behind  such  a  detonation  wave  is  constant 
along  lines  of  r/t  =  constant.  One  can  then  define ‘similarity 
variables 


X  =  r/r^  or  T  =  t/t^  (2) 

and  the  flow  expands  self-sirailarly,  i.e.,,  it  is  constant 
along  lines  of  x  «  constant  or  t  ■  constant.  The  analysis 
of  self-similar  blast  wave  flow  fields  is  performed  most 
easily  for  the  so-called  phase  plane  variables: 

P  «  u/xw^  «  TU/W^  (3) 

n  n 

Z  »  (a/xw^)^  »  (xa/w  (4) 

n  H 

where  a  de.notes  the  isentropic  sound  speed,  a  >  /yp/i  .  Por 
the  above  assumptions,  the  conservation  laws  reduce  to  a 
single  ordinary  differential  equation  (Kuhl  and  Seirew,  1978) s 

(5) 

(6) 


^  3  20+j (Y-1) (l-PiP 

ar  •  f  — - 

D  • 


and  a  quadrature 
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The  thermodynamic  behavior  of  the  medium  is  characterized  by 
the  specific  heat  ratio,  y,  which  is  assumed  to  be  constant. 
Flow  geometry  is  taken  into  account  by  the  paramter  j  which 
equals  0,  1  or  2  for  plane,  line  or  point  symmetric  flow. 

Equations  (5)  and  (7)  are  integrated  from  the  wave  front 

(F  ,  z  ,  X  =  T  =  1)  to  the  center  of  the  wave.  Given  the 
n  n 

integral  curve  Z  =  Z(F)  and  x  =  x(P)  (or  x  =  t(F)),  the  flow 
field  is  determined  from 


u/u^  =  xF/F^ 

(8) 

=  F/xF^ 

e/e^  =  (p/Pj^)  =  (p/pj^)^“^  =  ^^2/^n  “ 

the  latter  being  a  consequence  of  the  fact  that  the  flow 
field  behind  the  detonation  wave  is  homentropic. 


2 .  BOUNDARY  CONDITIONS 


The  physical  variables  just  behind  the  detonation  front 
are  determined  from  conservation  of  mass,  momentum  and  energy 
across  the  front: 


Pn  “  Pa  * 

2  2 

n  a  n  n 


(10) 
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\vritten  here  in  terms  of  F^,  and  the  chemical  energy  release, 
q.  The  CJ  state,  where  the  wave  velocity  'equals  the  speed 
of  the  characteristic 


w 


n 


(11) 


then  becomes 


n  n 

=  (l-y)/(Y+l)  (12) 

2  2 

where  y  is  the  front  Mach  number  parameter  y  = 

Typical  detonation  Mach  numbers  are  on  the  order  of  M  «  5 
(y  »  0.04)?  hence,  the  strong  shocJc  boundary  conditions 
(y  =  0)  are  a  good  approximation  for  most  problems.  The 
phase  plane  variables  at  the  detonation  front  are  then 


1/(Y+1) 


-  Iy/(y+1)]  , 


while  the  jump  conditions  reduce  to 


®n  * 


p  ■  p  +p  /(y+1) 


«  -  1 
'’n  ^ 
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(13) 


(14) 


(14) 


k  =  /2(Y^'l)  . 

As  Shown  above,  all  properties  behind  strong 
as  well  as  the  heat  release  at  the  front  are 
the  detonation  velocity  w^  and  y. 

3.  PLANAR  DETONATIONS 

Planar  detonations  (j  ®  0)  represent  a  singular  case  in 
which  the  solution  may  be  written  in  closed  form.  The  rare¬ 
faction  wave  behind  the  detonation  front  is  a  centered, 
simple  wave  which  decelerates  the  flow  to  zero  velocity 
(u  =  0  at  r  =  0).  One  finds  from  characteristics  theory  that 
along  lines  of 

dr 

■  w  ■  u+a 

•  xw„[P+2’®),  (15) 

the  Riemann  variable 

R  »  m*2a/(y-l) 

*  xw„lF-22V(Y-l)]  (16) 

n 

remains  constant.  Equation  (IS)  can  be  solved  (using  the 
fact  that  w  »  xw^)  to  obtain  the  integral  curve  in  the  phase 
plane  in  closed  form: 


CJ  detonations 
functions  of 


2  •  (1-P)^. 
11 
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Eliminating  Z  from  Equation  (16)  by  virtue  of  Equation  (17), 
one  find!? 


R  =  xw^[(Y+1)F-2]/(y-1).  (18) 

Evaluating  the  Riemann  constant  of  Equation  (18)  at  the 
detonation  front,  yields 

xw^[(Y+l)F-2l/(Y-l)  =  W^[(Y+l)F^-2]/(Y-l);  (19) 

and  solving  for  x  one  finds 


2-(Y+l)Fp 

^  2-(Y+1)F  ' 

(20) 

or,  by  virtue  of  Equation  (7),  x  becomes; 

2-(y+1)F 

T  “  Vx  -  2=  (y+1T^„  • 

(21) 

These  re{..esent  an  analytic,  closed-form 
integral  of  Equation  (7).  For  the  strong 
condition  ( (•''^  =  1/(Y+1)},  Equation  (21) 
particularly  simple  fCi.m: 

relation  for  the 

shoc}c  boundary 
acquires  the 

T  »  2-F{y+1). 

(22) 

Equations  (17)  and  (22)  can  oe  used  to  express  the  phase 
plane  variables  in  terms  of  x: 

F  »  {2-x)/{y+1) 

(23) 

2  *  {y+x-1)  V:y+1)^. 

(24) 

12 


i 


I 


The  above  relations  combined  with  the  boundary  conditions 
[Equations  (13)  and  (14)]  can  be  used  in  Equations  (8)  and' 

(9)  to  give  an  analytic  expression  for  the  flow  field  of  a 
strong  planar  CJ  detonation  wave: 


u/Un  = 

2/t-I 

(25) 

a/.n  = 

[ (y+t-1)/yt]^ 

(26) 

P/P„  = 

[(Y+t-D/yt]^^''^^"^^ 

(27) 

P/P„  - 

t(Y+T-l)/YT]^/^^"^^ 

(28) 

q/q„  . 

(2/t-I)  ^  [  (Y+T-1)  /Yt] 

(29) 

This  solution  is  valid  from  the  head  to  the  tail  of  the 
rarefaction  wave:  1  <  t  <  2.  For  times  greater  than  t  =  2, 
the  flow  field  is  constant: 


u/u^  •  0 
n 

(30) 

e/Sn  •  (Y+1)V4y^ 

(31) 

The  above  restrictions  allow  the  solution  to  satisfy  the 
u  ••  0  condition  at  the  center  of  symmetry,  x  •  0. 

Figure  3  shows  the  derived  similarity  solution  for  a 
constant-velocity  CJ  detonation  wave  (j  «  0,  y  *  0,  y  *  1.3) 
as  a  function  of  time  t,  with  all  quantities  being  nondimen- 
sionalized  by  their  values  at  the  detonation  front.  An 
auxiliary  x  scale  is  included  to  display  the  space  distri¬ 
bution  aspects  of  the  flow  field.  Figure  4  gives  the 
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Figure  3.  Flow  field  time  history  for  a  planar 

Chapman-Jouguet  detonation  wave  (j  »  0, 
Y  “  1*3»  y  *  0) . 


Figure  4 .  Nondimens ional  static  pressure  impulse  and 
dynamic  pressure  impulse  Z2  ^ 

planar  detonation  wave  (j^0^Y*l*3#y»0). 
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* 

nondimens ional  static  and  dynamic  pressure  impulses 


13^(1)  = 


32) 


l2(T) 


(33) 


as  a  function  of  time  t.  At  the  end  of  the  “positive  phase" 

(t  =  2),  the  integrals  reach  the  following  values  for  the 
case  of  Y  =  1.3:  =  2)  =  0.5513  and  I2(t  *  2)  =  0.1818. 

This  figure  can  be  used  to  evaluate  the  overpressure  and 
dynamic  pressure  partial  impulses  from  the  following  relations: 


‘n 

y,  ^ 

Iq(t/tj^)  »  J  qdt 


Dynamic  pressure  defined  as  q 


2 

»  Hpu  . 
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III.  APPLICATIONS 


1.  PROPERTIES  OF  DETONATABLE  GASES 

Table  la  shows  the  detonation  properties  of  gaseous 
mixtures  as  originally  calculated  and  measured  by  Jouguet 
(1905-6).  With  improved  thermochemical  data,  Lewis  and 
von  Elbe  (1961)  calculated  the  detonation  properties  for 
hydrogen-oxygen  systems  given  in  Table  lb.  Table  Ic  presents 
detonation  velocities  for  hydro-carbon  plus  oxygen  or  air 
mixtures  from  Shchelkin  and  Troshin  (1965) .  Detonation 
velocities  for  hydrogen-oxygen  systems  range  from  1.7  to 
3.S  )an/s  depending  on  the  mixture,  while  detonation  pressures 
vary  from  *  14  to  18  atm^  detonation  temperatures  range 
from  2600  to  3600  K.  Acetylene  (CjH^)  or  cyanogen  (CjN^)  plus 
oxygen  mixtures  have  considerably  higher  detonation  pressures 
(^>55  bars)  and  temperatures  approaching  6000  K. 

2.  SIMULATION  FIDELITY 

As  a  representative  case,  we  consider  here  a  stoichiome¬ 
tric  mixture  of  hydrogen  and  oxygen  at  atmospheric 

pressure  and  temperature.  The  detonation  jump  conditions  are 
listed  in  Table  2.  Detonation  pressure  is  17.6  atm,  while 
the  peak  dynamic  pressure  and  temperature  are  6  atm  and 
about  3660  K,  respectively. 

Figure  S  compares  the  time  histories  for  a  planar,  CJ 
detonation  wave  solution  {j»0,  y*1.3,  y^O)  with  a 
nuclear  surface  burst  explosion  for  the  case  where  the  peak 
pressures  have  been  matched  at  p^^  ^  17.6  atm  »  244  psi) . 
Times  have  been  scaled  by  the  wave  arrival  time  at  the 
station  of  interest.  Notice  that  the  CJ  pressure  waveform 
decays  much  more  slowly  than  the  nuclear  case.  This  is  very 
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TABLE  2 


DETONATION  CONDITIONS  FOR  A 
STOICHIOMETRIC  H2-O2  MIXTURE 
(w  =  2.821  km/s) 


Detonation 

Air  shock 

Ambiant 

State  n 

1.0 

17.6 

17.8 

0.SO3  X  lO"® 

0.887  X  10"^ 

5.74  X  10*3 

291 

••3860 

291 

12 

^  15 

29 

1.4 

1.30 

1.37 

0.531 

1.824 

0.646 

0 

1.813 

1.04 

0 

6.067 

31.0 

- 

2J2 

1.32 
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Figure  S.  Comparison  of  rhe  firoe  histories  for  the  planar  chapman -Jouguet 
detonation  wave  solution  (CJ)  and  a  nuclear  surface  burst 
explosion  (PEI  for  the  case  of  matched  peak  pressures  (p  =  17. 


advantageous  because  it  will  reduce  the  required  tube  length 
for  the  simulator.  The  time  scale  for  the  nuclear  case  can 
be  stretched  to  fit  the  CJ  solution  according  to  the  relation; 


T  =  ^  2.224 

XJ  'PE 


(36) 


which  is  denoted  by  the  circled  points  on  Figure  5.  Hence, 
to  good  approximation,  one  finds: 

P/Pjj  =  ((0.3/T(,j+1)/1.31^®''^ 

=  t(0.3/Tpg^-^^^+l)/1.3]^®/^  (37) 

It  is  clear  from  the  above  comparisons  that  the  detonation 
tube  simulator  does  a  good  job  at  replicating  the  pressure 
waveforms  of  a  nuclear  surface  burst  for  times  1  <  t  <2.  The 
next  question  is:  how  well  does  it  match  the  other  S^viron* 
ments?  Figure  S  compares  the  dynamic  pressure  histories  for 
the  CJ  and  PE  cases.  The  waveforms  have  the  same  shape;  how¬ 
ever,  the  air  shock  has  a  much  higher  peak  dynamic  pressure 
(31  atm)  than  the  CJ  detonation  (6  atm).  This  is  caused 
principally  by  the  low  initial  density  of  the  H2*’02  mixture 
at  ambient  pressure  and  temperature  (which  is,  in  turn,  re¬ 
lated  to  its  smaller  molecular  weight  of  12  versus  29  for  air) . 
This  peak  dynamic  pressure  can  be  increased  by  choosing 
reactants  with  large  molecular  weights.  One  of  the  best 
mixtures  in  this  regard  is  a  cyanogen-oxygen  mixture 
(02^2^202)  which  gives  a  peak  dynamic  pressure  of  34.6  atm. 

There  are  significant  differences  in  the  temperature 
histories  for  the  CJ  and  PE  eases,  as  demonstrated  in  Fig¬ 
ure  S.  At  the  17.6-atm  station,  the  air  is  shock-heated  to 
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a  temperature  of  about  1000  K  and  increases  to  about  3000  K 
during  the  simulation  time  (1  <  i  <  2),  while  the  CJ  temperature 
history  starts  at  the  detonation  temperature  ('^'3700  K  for 
the  H2-O2  system  considered)  and  decays  to  about  3000  K  at 
T  =  2,  At  these  gas  temperatures  small  dust  particles  can  melt 
and/or  vaporize  due  to  convective  heat  transfer  between  the 
detonation  products  and  the  dust  (e.g.,  Si02  has  a  melt  temp¬ 
erature  of  1900  K  and  a  vaporization  temperature  of  2500  K). 
Although  air  temperatures  are  too  small  to  cause  this  effect 
in  the  nuclear  case  (until  a  time  of  t  =  1.6  to  1,8),  the  dust 
particles  will  be  heated  by  the  nuclear  fireball  radiation, 
and  may  melt  and/or  vaporize  due  to  radiative  heat  transfer. 

To  resolve  this  issue,  detailed  heat  transfer  calculations 
should  be  performed  for  dust  particles  of  interest.  It  will 
probably  be  difficult  to  find  gaseous  mixtures  with  comous- 
tion  temperatures  much  below  2500  K,  so  if  heat  transfer  is 
a  problem,  then  it  will  place  limits  on  the  range  of  applica¬ 
bility  of  the  simulator.  Alternatively,  one  could  match  the 
detonation  temperature;  a  3500  K  air  shock  temperature  cor¬ 
responds  to  a  shock  overpressure  of  about  100  atm. 

Current  theories  of  boundary  layer  scouring  of  dust  from 
surfaces  attempt  to  relate  mass  injection  rates  to  the  wall 
shear  stress.  This  is,  in  turn,  related  to  the  dynamic  pres¬ 
sure  just  outside  the  boundary  layer.  Hence,  the  dynamic 
pressure  appears  to  be  the  most  important  flow  parameter. 

Figure  6  compares  the  dynamic  pressure  time  histories  for 
the  detonation  tube  simulator  (using  the  ^2^0^  systett.  shown 
in  Figure  S)  with  those  of  the  PE  case  at  matched  peak 
dynamic  pressures  (q^  «  6  atm).  It  is  clear  from  Figure  6  that 
the  CJ  waveform  is  a  good  approximation  to  the  PE  waveform 
over  most  of  the  positive  phase;  hence,  for  scaling  purposes 

tcj  -  (38) 
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Figure  6.  Comparison  of  the  time  histories  for  the  planar  Chapman-Jouguet 
detonation  wave  solution  (CJ)  and  a  nuclear  point  explosion  (PE) 
for  the  case  of  matched  peak  dynamic  pressures  (q  =  6  atm) . 


for  dynamic  pressures.  The  analytic  form  for  the  dynamic 
pressure  history  then  becomes 

q/g^  =  (2/t-1)^[(0.3/t+1)/1.3]^°/^.  (39) 

In  this  case  the  peak  overpressures  are  no  longer  matched 

(Ap  =  16.6  atm  in  the  CJ  case,  and  Ap  5  atm  in  the  PE 
n  n 

case) .  Temperatures  are  also  mismatched,  as  before. 

3.  SCALING 


Next  let  us  consider  scaling  relations  for  the  detona¬ 
tion  tube  blast  wave  simulator.  We  start  by  equating  simu¬ 
lation  times; 


At 


CJ 


At. 


PE 


(40) 


or 


(41) 


Solving  for  t^^^^  and  noting  that  =  ip^  (where  a  ■  2.224 
for  static  pressure  wave  forms  and  a  ■■  1  for  dynamic  pressure 
wave  forms  according  to  Equations  (36)  and  (38) #  respectively), 
one  finds 


t  «  t  W^/^(T„--1)/(T„,,®‘-1) 
"CJ  "KT  PE  PE  ^ 


(42) 


where  «  shock  arrival  time  in  ms/KTgg^^^,  and 


W  »  equivalent  nuclear  burst  yield  in  kilotons. 
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The  above  relation  can  be  used  to  determine  the  detonation 
run  length  (L  =  needed  to  provide  the  required  simu¬ 

lation  time; 


L 


CJ  nKT 


-1). 


(43J 


This  distance  is  considered  to  be  the  appropriate  location 

of  the  dusty  flow  measure  devices.  As  shown  in  Table  3/  a 

1/3 

run  length  of  26  m/KTgg  is  required  to  simulate  the 
pressure  wave  form  at  the  p^^  =  17.6  atm  (Ap  =  244  psig)  sta¬ 
tion  for  a  kiloton  surface  burst  with  a  stoichiometric  mix¬ 
ture  of  hydrogen  and  oxygen,  with  this  same  mixture,  a  run 

1/3 

length  of  182  ra/KTgg  is  required  to  simulate  a  kiloton 
dynamic  pressure  wave  form  at  the  =  6  atm  (88  psi)  station. 
As  one  goes  to  smaller  pressures,  the  nuclear  waveforms 
broaden;  hence,  longer  run  lengths  are  required.  To  achieve 
larger  dynamic  pressures  one  can  use,  for  example,  an  acety¬ 
lene-oxygen  mixture  (C,H,+0-)  with  a  peak  dynamic  pressure  of 

^  ^  ^  1/3 

18  atm  and  a  run  length  of  lOS  m/KTgg  '  ,  or  a  cyanogen- 
oxygen  mixture  (C,N,+20,)  with  a  peak  dynamic  pressure  of 

^  ^  1/3 

35  atm  and  a  run  length  of  66  m/KTgg  '  . 

4.  DESIGN  CONSIDERATIONS 


Next  we  offer  the  following  practical  considerations  for 
fielding  a  blast  wave  simulator:  the  choce  of  gaseous  mix¬ 
ture,  a  technique  for  filling  the  detonation  tube,  and 
approaches  to  minimize  the  tube  length  (hence,  minimize  cost). 
The  current  costs  of  a  few  candidate  gases  are: 

gas:  S5  per  100  ft^  (at  2000  psi  and  291  K) 

Oj  gas:  $1.87  per  100  ft^  tat  2000  psi  and  291  K) 

CjHj  gas:  $14.75  per  100  ft^  (at  250  psi  and  291  K) 

C2N2  gas:  $340  per  lb  (2-in  dia  x  13-in  cylinder) 

2S 


Omous  limutaiion  duration 


-w 


Cyanogen  (C2N2)  appears  to  be  prohibitively  expensive. 

Assuming  a  1-m^  tube  cross  section,  it  would  cost  about  $5 
to  fill  a  115-ra-long  tube  with  acetylene  for  an  18-atm 
dynamic  pressure  simulator,  and  less  than  $1  to  fill  a  35-m- 
long  tube  with  2H2+O2  as  an  18-atm  static  pressure  simulator. 
Clearly  these  are  negligible  costs. 

Perhaps  the  simplest  way  of  filling  the  tube  with  a 
detonable  gas  mixture  is  to  use  air  as  the  oxidizer  and 
slowly  bleed  in  the  fuel  gas.  Utilizing  air  instead  of  O2 
gas  will  reduce  the  detonation  velocity  to  about  1.9  km/s 
(which  will  reduce  the  run  length  required  by  about  30  per¬ 
cent)  and  reduce  the  detonation  temperature  to  per.haps 
2700  K— both  of  which  are  beneficial  effects.  Also,  one 
would  like  to  choose  a  gas  mixture  which  has  a  high  molecular 
weight,  to  achieve  large  peak  dynamic  pressures.  Taking  all 
these  factors  into  consideration,  it  appears  that  a  stoichio¬ 
metric  mixture  of  acetylene  and  air  is  one  of  the  best  choices 
for  the  detonatable  gas.  It  is  inexpensive  ('v$5  per  100-m^ 
test),  it  gives  large  peak  dynamic  pressures,  and  it  is 
readily  available. 

Probably  the  largest  expense  item  of  the  simulator  is 
the  construction  cost  of  building  the  tube.  This,  in  turn, 
is  related  to  the  tube  length.  Simulation  time  requirements 
fix  the  detonation  run  length,  but  a  downstream  tube  section 
is  also  required.  This  end  of  the  tube  could  be  sealed  with 
a  thin  plastic  sheet,  and  terminated  with  a  survivable 
louvered  end  section  (Crosnier  et  al.,  1974}  to  eliminate 
rarefactions  from  this  "semi-open**  end  and  minimize  the  tube 
length. 
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IV.  CONCLUSIONS 


When  a  planar  detonation  wave  is  initiated  at  the  closed 
end  of  a  tube,  it  induces  a  peaked  blast  wave  which  expands 
self-similar ly  with  time.  Being  forced  by  the  physics,  this 
happens  naturally.  The  preceding  analysis  demonstrated  that 
the  static  and  dynamic  pressure  waveforms  associated  with  a 
planar  CJ  detonation  give  a  high-fidelity  simulation  of  either 
the  static  or  dyneunic  pressure  environment  (but  not  both  at 
the  same  time)  of  a  nuclear  surface  burst  for  times  1  <  <  2 

One  can  achieve  peak  overpressures  from  14  to  55  atm  and  peak 
dynamic  pressures  from  6  to  35  atm,  depending  upon  the  detona- 
table  gas  selected.  This  phenomenon  could  be  utilized  in  a 
disposable  detonation  tube  to  impose  peaked  blast  waves  on  in- 
situ  ground  surfaces  of  interest,  and  then  measure  the  proper¬ 
ties  of  the  ensuing  dusty  (or  clean)  blast  wave  boundary  layer. 
Detonation  tube  lengths  from  25  to  100  m/KTgg  ' *  are  required 
to  simulate  yields  in  the  kiloton  range. 

The  principal  deficiency  of  this  technique  is  that  the 
detonation  products  have  temperatures  ('v3000  K)  which  are 
significantly  larger  than  their  airbiast  counterparts  ('V'lOOO  K) 
Small  (10-iim  diameter)  dust  particles  have  the  potential  for 
melting  and  vaporizing  due  to  convective  heat  transfer.  These 
effects  can  be  ameliorated  somewhat  by  choosing  air  as  the 
oxidizer  and  by  using  large  diameter  dust  particles. 
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